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Symmetries play a crucial role in electroweak symmetry breaking models with non- 



minimal Higgs content. Within each class of these models, it is desirable to know which 



symmetry groups can be implemented via the scalar sector. In A-Higgs-doublet mod- 
els, this classification problem was solved only for N = 2 doublets. Very recently, we 
suggested a method to classify all realizable finite symmetry groups of Higgs-family trans- 
formations in the three-Higgs-doublet model (3HDM). Here, we present this classification 
in all detail together with an introduction to the theory of solvable groups, which play 
the key role in our derivation. We also consider generalized- CP symmetries, and dis- 
cuss the interplay between Higgs-family symmetries and CP-conservation. In particular, 
we prove that presence of the Z4 symmetry guarantees the explicit CP-conservation of 
the potential. This work completes classification of finite reparametrization symmetry 
groups in 3HDM. 



1 Introduction 



X 1 

The nature of the electroweak symmetry breaking is one of the main puzzles in high-energy 
physics. Very recently, the CMS and ATLAS collaborations at the LHC announced the dis- 
covery of the Higgs-like resonance at 126 GeV, [lj, and their first measurements indicate 
intriguing deviations from the Standard Model (SM) expectations. Whether these data signal 
that a non-minimal Higgs mechanism is indeed at work and if so what it is, are among the 
hottest questions in particle physics these days. 

In the past decades, many non- minimal Higgs sectors have been considered, [2J. One 
conceptually simple and phenomenologically attractive class of models involves several Higgs 
doublets with identical quantum numbers (NHDM, A^-Higgs-doublet models). Its simplest 
version with only two doublets, 2HDM, was proposed decades ago, [3], but it is still actively 
studied, see [I] for a recent review, and it has now become a standard reference model of the 
beyond Standard Model (bSM) physics. Constructions with more than two doublets are also 
extensively investigated, [5T4T4"]. 

Many bSM models aim at providing a natural explanation for the numerical values of 
(some of) the SM parameters, and it is often done by invoking additional symmetries in the 
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model. These are not related with the gauge symmetries of the SM but rather reflect extra 
symmetry structures in the "horizontal space" of the model. One of the main phenomenological 
motivations in working with several doublets is the ease with which one can introduce various 
symmetry groups. Indeed, Higgs fields with identical quantum numbers can mix, and it is 
possible that some of these Higgs-family transformations leave the scalar sector invariant. Even 
in 2HDM, presence of such a symmetry in the lagrangian and its possible spontaneous violation 
can lead to a number of remarkable phenomena such as various forms of CP- violation, [3|[T5]. 
non-standard thermal phase transitions which may be relevant for the early Universe, [T6] , 
natural scalar dark matter candidates, [17J. For models with three or more doublets, an extra 
motivation is the possibility to incorporate into the Higgs sector non-abelian finite symmetry 
groups, which can then lead to interesting patterns in the fermionic mass matrices (for a 
general introduction into discrete symmetry groups relevant for particle physics, see [H]). In 
this respect, the very popular symmetry group has been A4, [HIE], the smallest finite group 
with a three-dimensional irreducible representation, but larger symmetry groups also received 
some attention, [5}l9j[TU]. 

Given the importance of symmetries for the NHDM phenomenology, it is natural to ask: 
which symmetry groups can be implemented in the scalar sector of NHDM for a given N? 

In the two-Higgs-doublet model (2HDM), this question has been answered several years 
ago, [19j[20], see also [I] for a review. Focusing on discrete symmetries, the only realizable 
Higgs-family symmetry groups are Z2 and (Z2) 2 . The Z2 group can be generated, for example, 
by the sign flip of one of the doublets (and it does not matter which, because once we focus 
on the scalar sector only, the simultaneous sign flip of both doublets does not change the 
lagrangian), while the (Z2) 2 group is generated by sign flips and the exchange <pi <H- <p2- If 
generalized- CP transformations are also included, then (Z2) 3 becomes realizable as well, the 
additional generator being simply the CP conjugation. 

With more than two doublets, the problem remains open. Although several attempts have 
been made in past to classify at least some symmetries in NHDM, [TTJESldl] > they led only 
to very partial results. The main obstacle here was the lack of the completeness criterion. 
Although many obvious symmetry groups could be immediately guessed, it was not clear how 
to prove that the given potential does not have other symmetries. An even more difficult 
problem is to prove that no other symmetry group can be implemented for a given N . 

In the recent paper [21] we found such a criterion for abelian symmetry groups in NHDM for 
arbitrary N. Since abelian subgroups are the basic building blocks of any group, classification 
of realizable abelian symmetry groups in NHDM was an important milestone. We stress that 
this task is different from just classifying all abelian subgroups of SU(3), because invariance 
of the Higgs potential places strong and non-trivial restrictions on possible symmetry groups. 

In this paper, we solve the classification problem for all finite symmetry groups in 3HDM, 
including non-abelian groups. We do this by using the abelian groups in 3HDM found in [21] 
and by applying certain results and methods from the theory of solvable groups. Some of these 
results were already briefly described in [22]. Here, we present a detailed derivation of this 
classification together with an introduction to the relevant methods from finite group theory. 
In addition, we extend the analysis to symmetry groups which include both Higgs-family and 
generalized- CP transformations. This work, therefore, solves the problem of classification of 
finite reparametrization symmetry groups in 3HDM. 

We would like to stress one important feature in which our method differs from more 
traditional approaches to symmetry classification problem, at least within the bSM physics. 
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Usually, one starts by imposing invariance under certain transformations, and then one tries 
to recognize the symmetry group of the resulting potential. In this way it is very difficult to 
see whether all possible symmetries are exhausted. We approach the problem the other way 
around. We first restrict the list of finite groups which can appear as symmetry groups of 
3HDM, and then we check one by one whether these groups can indeed be implemented. 

The structure of this paper is the following. In Section [2] we describe different types of 
symmetries in the scalar sector of NHDM and discuss the important concept of realizable 
symmetry groups. Section [3] contains an elementary introduction into the theory of (finite) 
solvable groups. Although it contains pure mathematics, we put it in the main text because 
it is a key part of the group-theoretic step of our classification, which is presented in Sec- 
tion HI Then, in Section [5] we describe the methods which we will use to prove the absence 
of continuous symmetries. Sections [6] and [7] contain the main results of the paper: explicit 
constructions of the realizable symmetry groups and of the potentials symmetric under each 
group. Finally, in Section |H] we summarize and discuss our results. For the reader's con- 
venience, examples of potentials for each of the realizable non-abelian symmetry groups are 
collected in the Appendix. 

2 Symmetries of the scalar sector of multi-Higgs-doublet 
models 

2.1 Reparametrization transformations 

In NHDM we introduce N complex Higgs doublets with the electroweak isospin Y = 1/2, 
which interact with the gauge bosons and matter fields in the standard way, and also self- 
interact via a Higgs potential. The generic renormalizable Higgs potential can contain only 
quadratic and quartic gauge-invariant terms, and it can be compactly written as [231121] : 

V = Y ab (^J b ) + Z abcd (<P%M<Pd) , (1) 

where all indices run from 1 to N. Coefficients of the potential are grouped into components 
of tensors Y ab and Z abcd ; there are N 2 independent components in Y and N 2 (N 2 + l)/2 
independent components in Z. 

In this work we focus only on the scalar sector of the NHDM. Therefore, once coefficients 
Y ab and Z abc d are given, the model is completely defined, and one should be able to express 
all its properties (the number and the positions of extrema, the spectrum and interactions of 
the physical Higgs bosons) via components of F's and Z's. This explicit expression, however, 
cannot be written via elementary functions, and it remains unknown in the general case for 
any N > 2. 

A very important feature of the most general potential is that any non-degenerate lin- 
ear transformation in the space of Higgs doublets belonging to the group GL(2, C) keeps the 
generic form of the potential, changing only the coefficients of Y and Z. We call such a trans- 
formation a Higgs-basis change. In addition, the CP transformation, which maps doublets 
to their hermitean conjugates <p a — > (f)' a , also keeps the generic form of the potential, up to 
coefficient modification. Its combination with a Higgs-basis change represents a transforma- 
tion which is usually called a generalized-CP transformation, [25]. The Higgs basis changes 
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and generalized- CP transformations can be called together reparametrization transformations 
because they preserve the generic structure of the potential and lead only to its reparametriza- 
tion. 

A reparametrization transformation changes the basis in the space of Higgs doublets but 
does not modify the structural features of the model such as the number and the properties 
of minima, the symmetries of the potential and their spontaneous breaking at the minimum 
point. These properties must be the same for all the potentials linked by reparametrization 
transformations. Therefore, these properties must be expressible in terms of reparametrization- 
invariant combinations of Y's and Z's, [231126] . 

If a reparametrization transformation maps a certain potential exactly to itself, that is, 
if it leaves certain F's and Z's invariant, we say that the potential has a reparametrization 
symmetry. Usually, there is a close relation between the reparametrization symmetry group G 
of the potential and its phenomenological properties, both within the scalar and the fermion 
sectors. Therefore, understanding which groups can appear as reparametrization symmetry 
groups in NHDM with given N is of much importance for phenomenology of the model. 

2.2 The group of kinetic-term-preserving reparametrization trans- 
formations 

Often, one restricts the group of reparametrization transformations only to those transfor- 
mations which keep the Higgs kinetic term invariant. In this case, a generic basis change 
becomes a unitary transformation <ft a y U^ipb with U G U(N). A kinetic-term-preserving 
generalized- CP transformation is an anti-unitary map <p a h-> U a b<f)l, which can be written as 
Ucp = U ■ J, with a unitary U and with J being the symbol for the CP-transformation. 

The group U(N) contains the subgroup of overall phase rotations, which is already included 
in the gauge group U(1)y- Since we want to study structural symmetries of the NHDM 
potentials, we should disregard transformations which leave all the potentals invariant by 
construction. This leads us to the group U(N)/U(1) ~ PSU(N). Note that SU(N), which 
is often considered in these circumstances, still contains transformations which only amount 
to the overall phase shift of all doublets. They form the center of SU(N), Z(SU(N)) ~ Z N , 
and act trivially on all NHDM potentials. Being invariant under them does not represent 
any structural property of the Higgs potential, therefore, we are led again to the factor group 
SU(N)/Z(SU(N)) = PSU(N). This allows us to write the group of kinetic-term-preserving 
reparametrization transformations as a semidirect product of the Higgs basis change group 
and the Z 2 group generated by J (for a more detailed discussion, see [21]): 

G rep = PSU(N) x Z* . (2) 

Here the asterisk indicates that the generator of the corresponding group is an anti-unitary 
transformation; we will use this notation throughout the paper. 

Below, when discussing symmetry groups of the 3HDM potential, we will be either looking 
for subgroups of PSU(3) (if only unitary transformations are allowed) or subgroups of this 
Grep (when anti- unitary reparametrization transformations are also included). This should 
always be kept in mind when comparing our results with the groups which are discussed 
as symmetry groups in the 3HDM scalar sector. For example, in [9],[T0] a 3HDM potential 
symmetric under A (27) or A(54) was considered, both groups being subgroups of S77(3). 



4 



However, they contains the center of SU(3), which, we repeat, acts trivially on all Higgs 
potentials. Therefore, the structural properties of that model are defined by the factor groups 
A(27)/Z(SU(3)) ~ Z 3 x Z 3 and A(54)/Z(SI/(3)) ~ (Z 3 x Z 3 ) x Z 2 , which belong to PSU(3). 

2.3 Realizable symmetry groups 

There is an important technical point which should be kept in mind when we classify symmetry 
groups of NHDM. When we impose a reparametrization symmetry group G on the potential, 
we restrict its coefficients in a certain way. It might happen then that the resulting potential 
becomes symmetic under a larger symmetry group G properly containing G. 

One drawback of this situation is that we do not have control over the true symmetry 
properties of the potential: if we construct a G-symmetric potential, we do not know a priori 
what is its full symmetry group G. This might be especially dangerous if G is finite while 
G turns out to be continuous, as it might lead to unwanted goldstone bosons. Another 
undesirable feature is related with symmetry breaking. Suppose that we impose invariance 
of the potential under group G but we do not check what is the true symmetry group G. 
After electroweak symmetry breaking, the symmetry group of the vacuum is G v < G, and it 
can happen that G v is not a subgroup of G. This is not what we normally expect when we 
construct a G-symmetric model, and it is an indication of a higher symmetry. 

Examples of these situations were encountered in literature before. For instance, the 
authors of [12] explicitly show that trying to impose a Z p , p > 2, group of rephasing transfor- 
mations in 2HDM unavoidably leads to a potential with continuous Peccei-Quinn symmetry. 
For 3HDM they find an even worse example, when a cyclic group immediately leads to a 
U(l) x L^l) -symmetric potential. Another well-known example is the v4 4 -symmetric 3HDM 
potential, which at certain values of parameters admits vacua with the S3 symmetry, although 
S3 is not a subgroup of A4, see an explicit study in jH]. The explanation is that the potential 
at these values of parameters becomes symmetric under S4 which contains both A4 and S3. 

In order to avoid such situations altogether, we must always check for each G whether 
the G-symmetric potentials are invariant under any larger group. We are interested only in 
those groups G, for which there exists a G-invariant potential with the property that no other 
reparametrization transformation leaves it invariant (either within PSU(3) or within G rep , 
depending on whether we include anti-unitary transformations). Following [14,21], we call 
such groups realizable. 

Using the terminology just introduced we can precisely formulate the two main questions 
which we address in this paper: 

1. considering only non-trivial kinetic-term-preserving Higgs-basis transformations (i.e. group 
PSU(3)), what are the realizable finite symmetry groups in 3HDM? 

2. more generally, considering non-trivial kinetic-term-preserving reparametrization trans- 
formations, which can now include generalized- CP transformations (i.e. group G rep ), 
what are the realizable finite symmetry groups in 3HDM? 

For abelian groups, these questions were answered in [21] for general N. Here we focus on 
non-abelian finite realizable groups for N = 3. 
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3 Solvable groups: an elementary introduction 



Our classification of realizable groups of Higgs-family symmetries in 3HDM contains two essen- 
tial parts: the group-theoretic and the calculational ones. The group-theoretic part will make 
use of some methods of pure finite group theory, which are not very familiar to the physics 
community (although they are quite elementary for a mathematician with expertise in group 
theory). To equip the reader with all the methods needed to understand the group-theoretic 
part of our analysis, we begin by giving a concise introduction to the theory of solvable groups. 
In doing so, we mention only methods and results which are relevant for the particular problem 
of this paper. For a deeper introduction to solvable groups and finite group theory in general, 
see e.g. [27] . 

3.1 Basics 

We assume that the reader is familiar with the basic definitions from group theory. We only 
stress here that we will work only with finite groups, therefore the order of the group G (the 
number of elements in G) denoted as |G| is always finite, and so is the order of any element g 
(the smallest positive integer n such that g n = e, the identity element of the group). 

A group G is called abelian if all its elements commute. An alternative way to formulate it 
is to say that all commutators in the group are trivial: [x, y] = xyx~ l y~ l = e for all x,y G G. 
Working with commutators is sometimes easier than checking the commutativity explicitly. 
For example, it is easy to prove that if every non-trivial element of the group has order two, 
g 2 = e, then the group is abelian. Indeed, for any x, y G G we have 

[x, y\ = xyx~ l y~ x = xyxy = (xyf = e , (3) 

which means that x and y commute. 

A group G can have proper subgroups H < G (whenever we do not require that the 
subgroup H is proper, we write H < G), whose order must, by Lagrange's theorem, divide the 
order of the group: \H\ divides \G\. If proper subgroups exist, some of them must be abelian. 
A simple way to obtain an abelian subgroup is to pick up an element g G G and consider its 
powers: if order of the element g is n, we will get the cyclic group Z n < G. 

The inverse of Largrange's theorem is not, generally speaking, true: namely, if p is a divisor 
of |G|, the group G does not necessarily have a subgroup of order p. However, if p is a prime 
which enters the prime decomposition of \G\, then according to Cauchy's theorem such a 
subgroup must exist (this group is Z p because there are no other groups of prime order). It 
immediately follows that if we have the list of all abelian subgroups of a given finite group G, 
then the prime decomposition of \G\ can only contain primes which are present in the orders 
of these abelian subgroups. 

In fact, there is an existence criterion stronger than Cauchy's theorem. Namely, if p a is 
the highest power of the prime p that enters the prime decomposition of \G\, then G contains 
a subgroup of this order, which is called the Sylow p-subgroup of the group G. This theorem 
(known as the Sylow-E theorem) is the starting point of the theory of Sylow subgroups, see 
Chapter 1 in [27|. 

There are several ways to present a finite group. One possibility is to list all its elements and 
write down the \G\ x \G\ multiplication table. Clearly, this presentation becomes impractical 
for a sufficiently large group. A more compact and powerful way is known as presentation by 
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generators and relations. We call a subset M = {gi, g 2 , . . . } of the elements of G a generating 
set (and its elements are called generators) if every g G G can be written as a product of 
elements of M or their inverses. The fact that G is generated by the set M is denoted as 
G = (M). Finding a minimal generating set for a given group and listing equalities which 
these generators satisfy is precisely presentation of the group by generators and relations. For 
example, the symmetry group of the regular n-sided polygon has the following presentation 
by generators and relations: 

D 2n = (a, b | a 2 = b 2 = (ab) n = e) . (4) 

This group is known as the dihedral group and has order \D 2n \ = 2n (note that there exists an 
alternative convention for denoting dihedral groups: D n ; the one which we use has its order 
in the subscript). 

3.2 Normal subgroups and extensions 

Consider two groups G and H. Suppose we have a map / from G to H, f : G — > H, which 
sends every g G G into its image f(g) G H. If this map preserves the group operation, 
f {9i)f{92) = f {9192)1 then it is called a homomorphism. If this map is surjective (i.e. it covers 
the entire H) and injective (distinct elements from G have distinct images in H), then / is 
invertible and is called an isomorphism. 

In the case when H = G, we deal with an isomorphism of the group onto itself, which 
is called an automorphism. One can note that composition of two automorphisms is also 
an automorphism, and define the group structure on the set of all automorphisms of G. 
This automorphism group is denoted as Aut(G). The trivial automorphism which fixes every 
element of G is the identity element of Aut(G). 

Let us now consider a special class of automorphisms called inner automorphisms, or 
conjugations. Fix an element g G G and define / : x !->■ g~ x xg for every x G G. It can 
be immediately checked that / is an automorphism, and that it sends a subgroup of G into 
a (possibly another) subgroup of G. It can however happen that certain subgroups will be 
mapped onto themselves: g~ l Hg = H . Subgroups which satisfy this invariance criterion for 
every possible g G G are called normal, or invariant subgroups. The fact that H is a normal 
subgroup of G is denoted as H <\G. 

Even when a subgroup H is not normal in G, one can pick up some elements g G G such 
that g~ 1 Hg = H . The set of elements of G with the property g~ l Hg = H forms a group, which 
is called the normalizer of H in G and denoted as N G (H). We then have H < N G (H) < G. 
Working with normalizers is a useful intermediate step in situations when it is not known 
whether the subgroup H is normal in G. 

Having a normal subgroup H<\G gives some information about the structure of G. One can 
define the group structure on the set of (left) cosets of H, which is now called the factor group 
G/H. Thus, one breaks the group into two smaller groups, which often simplifies its study. 
Given a normal subgroup H < G, one can define the canonical homomorphism : G — > G/H 
which sends every element g G G into its coset gH. Its kernel (all elements g which are 
mapped by into the identity element of G/H) is precisely H. Thus, every normal subgroup 
is the kernel of the corresponding canonical homomorphism. The reverse statement is also 
true: kernels of homomorphisms are always normal subgroups. 
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The group-constructing procedure inverse to factoring is called extension. Given two 
groups, N and H, a group G is called an extension of H by N (denoted as N .H), if there 
exists N < G such that N ~ N and G/N ~ H . In the case when, in addition, H is also 
isomorphic to a subgroup of G and G = NH, we deal with a split extension. The criterion 
for G to be a split extension can also be written as existence of N < G and H < G such that 
NH = G and N C\H — 1, so that G/N = H. The group G is then called a semidirect product 
G = N y\ H . 

Even if two groups N and if are fixed, they can support several extensions and split 
extensions. Therefore one faces the problem of classifying of all extensions of two given groups. 

For the most elementary example, consider extensions of H — Z 2 (generated by a) by 
iV = Z 2 (generated by b), which should produce a group of order 4. Then, for a split extension, 
we need a group G which has two distinct subgroups isomorphic to N and H. The only choice 
is G = Z 2 x Z 2 , which can be presented as (a, b | a 2 = b 2 = (ab) 2 = e). For a non-split 
extension, we require that only N is isomorphic to a subgroup of G. Thus, we still have 
b 2 = e, while a 2 must not be the unit element. Then we have to set a 2 = b producing the 
group Z 4 . So, Z 4 does not split over Z 2 , while Z 2 x Z 2 does. 

3.3 Characteristic subgroups 

In what concerns embedding of groups, normality is a relatively weak property. Namely, if 
K <\H and H <G, then K is not necessarily normal in G (it is instead called subnormal in G). 
Indeed, recall that a normal subgroup K < H stays invariant under all inner automorphisms 
on H. Here "inner" is meant with respect to the group H, namely, h~ 1 Kh = K for all 
h G H. However since H <\G, one can fix g G G but g H and consider an automorphism 
on H defined by H — > g~ l Hg. This is indeed an automorphism on H because it induces a 
permutation of elements of H preserving its group property, but it is not inner, because g does 
not belong to H . Therefore K does not have to be invariant under it: g~ x Kg 7^ K. 

However there is a stronger property which guarantees normality for embedded groups. 
Let us call a subgroup K characteristic in H if it is invariant under all (not only inner) 
automorphisms of H . Then, repeating the above arguments, we see that if K is characteristic 
in H, and H is normal in G, then K is also normal in G. Also, if K is characteristic in H 
and H is characteristic in G, then K is also characteristic in G. Thus, knowing that some 
subgroups are characteristic gives even more information than their normality. 

There is one simple rule which guarantees that certain subgroups are characteristic. If we 
have a rule defined in terms of the group G which identifies its subgroup H uniquely, then H 
is characteristic in G. Two important examples are: 

• the center of the group G denoted as Z(G), which is the set of all elements z G G such 
that they commute with all elements of G: 



The center of an abelian group coincides with the group itself. 

• the commutator subgroup (or derived subgroup) of G denoted as G' and defined as the 
subgroup generated by all commutators: 



Z(G) = {zeG\iz,g] = eVgeG}. 



(5) 



G' = ([x,y]), x,yeG. 



(6) 
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Note that the word "generated" is needed because the set of commutators is generally 
speaking not closed under the group multiplication. Clearly, the commutator subgroup 
of an abelian group is trivial, therefore the size of G can be used to qualitatively char- 
acterize how far G is from being abelian. 

3.4 Consequences of existence of a normal maximal abelian sub- 
group 

Let us now prove a rather simple group-theoretic result, which however will be important for 
our classification of symmetries in 3HDM. This result, loosely speaking, is the observation 
that mere existence of a subgroup of G with some special properties can strongly restrict the 
structure of the group G. 

First, an abelian subgroup A < G is called a maximal abelian subgroup if there is no other 
abelian subgroup B with property A < B < G. Note that the word "maximal" refers not 
to the size but to containment. This definition does not specify a unique subgroup; in fact 
a group can have several maximal abelian subgroups. They correspond to terminal points in 
the partially-ordered tree of abelian subgroups of G. 

Suppose that A is an abelian subgroup of a finite group G. Elemens of A, of course, 
commute among themselves. But it can also happen that there exist other elements g G G, 
g A, which also commute with all elements of A. The set of all such elements is called the 
centralizer of A in G: 

C G (A) = {geG\lg,a] = eVaeA}. (7) 

It is easy to check that Cg{A) is a subgroup of G, and it can be non-abelian. The name 
"centralizer" refers to the fact that although A is not the center in G, it is the center in 
C G (A). 

Clearly, A < C G (A). If A is a proper subgroup of Cg(A), then it means that A is not 
a maximal abelian subgroup. Indeed, we take an element g G C G (A), g A, and consider 
another subgroup B = (A,g). This subgroup is abelian and is strictly larger than A: A < 
B < G. On the other hand, an element x G G which commutes with all elements of B will 
certainly commute with all elements of A, while the inverse is not necessarily true. Therefore, 
we get the following chain: A < B < C G (B) < C G (A). Next, we check whether B is a proper 
subgroup of Cq{B). If so, we can enlarge it again in the same way by considering C = (B, g'), 
where g' G Cg(B), g' G" B. We can continue this procedure until it terminates with an abelian 
subgroup K which is self- centralizing: 

A < B < ■ ■ ■ < K = C G {K) <■■■< C G (B) < C G (A) . (8) 

Since there exists no other element in G which would commute with all elements of K, we 
conclude that K is a maximal abelian subgroup in G. 

Let us now see what changes if the abelian subgroup A is normal. Any element g G G 
acting on A by conjugation induces an automorphism of A. Thus, we have a map from G to 
the group of automorphisms of A, f : G — > Aut(A). The kernel of / consists of such g's which 
induce the trivial automorphism of A, that is, which leave every a G A unchanged: g~ 1 ag = a 
Va G A. But this coincides with the definition of centralizer. Therefore we conclude that 
kerf = C G (A). 
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The fact that Cq(A) is the kernel of the homomorphism / implies that Ca(A) is a normal 
subgroup of G. Note that it is essential that the abelian subgroup in question, A, is normal; 
if it were not, Cq{A) would not have to be normal. 

Now, if A is a normal maximal abelian subgroup of G, then ker/ = Cg(A) = A. In other 
words, the kernel of G/A — > Aut(A) is trivial, and therefore, G/A is isomorphic to a subgroup 
of Aut(A). Summarizing our discussion, if A is a normal maximal abelian subgroup of G, then 
G can be constructed as an extension of A by a subgroup of Aut(A): 

G-A.K, where K < Aut(A) . (9) 

This is a powerful structural implication for the group G of existence of a normal maximal 
abelian subgroup. 

3.5 Automorphism groups 

For future reference, we give some details on the automorphism groups Aut(A) of certain 
abelian groups A. In this subsection we will use the additive notation for the group operation. 
Suppose A = Z n is the cyclic group of order n with generator e: ne = e + ■ ■ ■ + e = 0. 

n times 

An automorphism a acting on A is a group-structure-preserving permutation of elements of 
A. Since A is generated by e, this automorphism is completely and uniquely defined once we 
assign the value of a(e) = k and make sure that ma(e) ^ for all < m < n. This holds 
when k and n are coprime (k — 1 is coprime to any n). The number of integers less than n 
and coprime to n is called the Euler function (p(n). Thus, we have \Aut{Z n )\ = (p(n). For a 
prime p, the Euler function is obviously ip(p) = p — 1. In general, if p^ 1 ■ ■ -p k s s is the prime 
decomposition for n, then 

<p(£ ■■■p k /) = M 1 ) ■ ■ ■ <f{p k s s ) = bi 1 - pi"- 1 ) • • • (p k s s - p'r 1 ) ■ 

Suppose now that p is prime and 

A = Z p x • • • x Z p = (Z p ) n . 

v ' 

n times 

Then G can be considered as an n-dimensional vector space over a finite field F p of order p. 
Vectors in this space can be written as 

x = k\t\ + • • • + k n e n , 

where numbers fcj G ¥ p and "basis vectors" are certain non-zero elements of the i-th group 
Z p . The group of all automorphisms on (Z p ) n is then the general linear group in this space 
GL n (p). 

Again, in order to define an automorphism a acting on A, it is sufficient to assign where 
the basis vectors are sent by a and to make sure that they stay linearly independent: that 
is, if m\o(e\) + • • • + m n a(e n ) = 0, with to* G F p , then all m,i = 0. In order to calculate 
\GL n (p)\, we just need to find to how many different bases the initial basis {ei, . . . , e n } can 
be mapped to. The first vector, ei, can be sent to p n — 1 vectors, the second vector, e%, can 
be then sent to p n — p vectors linearly independent with a(e 1 ), and so forth. The result is 

\GL n {p)\ = (p n - l){p n - p) ■ ■ ■ (p n - p n - 1 ) = p^(p - l)(p 2 - 1) • • • (p n - 1) . (10) 
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In particular, \Aut(Z p xZ p )\ = \GL 2 (p)\ = p(p — 1)Q» 2 — 1), and the p-subgroup of Aut{Z p x Z p ) 
can only be Z p . 

3.6 Nilpotent groups 

In group theory, a powerful tool to investigate structure and properties of groups is to establish 
existence of subgroup series with certain properties. For example, a finite collection of normal 
subgroups Ni < G is called a normal series for G if 

1 = iV < Nt < N 2 < ■ ■ ■ < N r = G . (11) 

Restricting the properties of the factor groups A^/iVj_x for all i, one can infer non-trivial 
consequences for the group G. 

If all the factor groups in the normal series lie in the centers, Ni/N^i < Z(G/Ni_i) for 
1 < % < r, then (ITT]) becomes the central series, and the group G is then called nilpotent. The 
smallest number r for which the central series exists is called the nilpotency class of G. 

Clearly, abelian groups are nilpotent groups of class 1 because for them G < Z(G). A 
non-abelian group G whose factor group by its center G/Z{G) gives an abelian group is a 
nilpotent group of class 2, etc. So, nilpotent groups are often regarded as "close relatives" of 
abelian groups in the class of non-abelian ones. One important class of nilpotent groups is 
p-groups, i.e. finite groups whose order is a power of a prime p. 

Nilpotent groups bear several remarkable features. We mention here only two of them 
which we will use below. First, a nilpotent group has a normal self-centralizing abelian sub- 
group (Lemma 4.16 in [27]), whose implications were discussed above. Second, if H is a proper 
subgroup of a nilpotent group G, then H is also a proper subgroup of N G (H) (Theorem 1.22 
in [27]). In other words, the only subgroup of a nilpotent group G which happens to be 
self-normalizing is the group G itself. 

3.7 Solvable groups 

A group G is called solvable if it has a normal series (ITT]) in which all factor groups Ni/N^i 
are abelian. This is a broader definition than the one of nilpotent groups. Therefore we can 
expect that both criteria and properties of solvable groups will be weaker than for nilpotent 
groups. 

One particular example is that unlike nilpotent grops, a solvable group does not have to 
possess a normal self-centralizing abelian subgroup. However what is does possess is just a 
normal abelian subgroup. In order to prove this statement, let us first introduce another series 
of nested subgroups, called the derived series. We first find G', the derived subgroup of G, 
then we find its derived subgroup, G" = (G 1 )', then the third derived subgroup, G^ = (G")', 
and so on. The derived series is simply 

■■• < G (3) < G" < G' < G. (12) 

The relation of the derived series with solvability is the following: G is solvable if and only 
if its derived series terminates, i.e. G^ = 1 for some integer m > (Lemma 3.9 in [27]). The 
basic idea behind the proof of this statement is the observations that G' is the unique smallest 
normal subgroup of G with an abelian factor group. Indeed, if iV < G and : G — > G/N is the 
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canonical homomorphism, then 4>{G') = {G/N)' (commutators are mapped into commutators). 
If we want G/N to be abelian, then (G/N)' = 1, and G' < ker</> = N. Therefore, whatever 
N r _i we choose in ffTTj) . it will contain G'. This argument can be continued through the series, 
and since the normal series terminates, so does the derived series. 
Now, since 

G (m) = 1 for 

some finite m, we can consider G (m_1) . It is an abelian group be- 
cause its derived subgroup is trivial. Being a characteristic subgroup of G ( - m ~ 2 \ it is definitely 
normal in G. Thus, we obtain the desired normal abelian subgroup. 

A normal abelian subgroup is not guaranteed to be maximal. One can, of course, extend it 
to a maximal abelian subgroup, but then it is not guaranteed to be normal. Thus, in order to 
use the result ([9]), we need to prove the existence of an abelian subgroup which combines both 
properties. This situation is not generic: a solvable groups does not have to possess a normal 
maximal abelian subgroup. However it can possess it in certain cases, and we will show below 
that in what concerns finite symmetry groups in 3HDM, they do contain such a group. 

4 Structure of the symmetry group in 3HDM 
4.1 Abelian subgroups and Burnside's theorem 

Our goal is to understand which finite groups G can be realized as Higgs-family symmetry 
groups in the scalar sector of 3HDM. We stress that we look for realizable groups only, see 
discussion in section 12.31 

Since finite groups have abelian subgroups, it is natural first to ask which abelian subgroups 
G can have. This can be immediately inferred from our paper [21] devoted to abelian symmetry 
groups in NHDM. In the particular case of 3HDM, only the following groups can appear as 
abelian subgroups of a finite realizable symmetry group G: 

Z 2 , Z 3 , Z 4 , Z 2 xZ 2 , Z 3 xZ 3 . (13) 

The first four are the only realizable finite subgroups of maximal tori in PSU(3). The last 
group, Z 3 x Z 3 , is on its own a maximal abelian subgroup of PSU(3), but it is not realizable 
because a Z3 x Z3-symmetric potential is automatically symmetric under (Z3 x Z3) x Z 2 , 
see explicit expressions below. However, since it appears as an abelian subgroup of a finite 
realizable group, it must be included into consideration. Trying to impose any other abelian 
Higgs-family symmetry group on the 3HDM potential unavoidably makes it symmetric under 
a continuous group. 

Let us first see what order the finite (non- abelian) group G can have. We note that the 
orders of all abelian groups in ( 1T3"j) have only two prime divisors: 2 and 3. Thus, by Cauchy's 
theorem, the order of the group G can also have only these two prime divisors: \G\ = 2 a 3 b . 
Then according to the Burnside's p a g b -theorem the group G is solvable (Theorem 7.8 in [27J), 
and this means that G contains a normal abelian subgroup, which belongs, of course, to the 
list (fl3l). 

In order to proceed further, we need to prove that one can in fact find a normal maximal 
(that is, self-centralizing) abelian subgroup of G, a property which is not generic to solvable 
groups but which holds in our case. 
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4.2 Existence of a normal abelian self-centralizing subgroup 

Suppose A < G is a normal abelian subgroup, whose existence follows from the solvability of 
G. In this subsection we prove that even if it is not self-centralizing, i.e. A < Cg{A), then 
there exists another abelian subgroup B > A, which is normal and self-centralizing in G. 

Suppose that A < Cg(A). Then for every b G Cg(A) \ A, the group Af, = (A,b) is an 
abelian subgroup of G, which properly contains A. There are two possibilities compatible with 
the list ([13]): 

(i) A = Z 2 , and then A^ can be either Z 2 x Z 2 or Z 4 , 

(ii) A = Z3, and then Af, = Z3 x Z3. 

Thus Cg{A) is either a 2-group or a 3-group. Below we assume that p — 2 if Cg{A) is a 
2-group, and p = 3 if Cg(^4) is a 3-group. 

Since Cg(^4) is a p-group, it is nilpotent, and according to discussion in section 13.61 it 
possesses a normal maximal abelian subgroup B, while B properly includes A = Z p : A < 
B < Cg(A). In particular, B is self-centralizing in Cg(A), so according to our discussion 
in section [331 the factor group Cq(A)/B is a subgroup of Aut(B). If B = Cg(A), then 
Cg(A) is abelian and, being a centralizer of a normal subgroup, it is normal in G. Clearly 
B < Cg{B) < Cg{A) = B, therefore Cg{A) is the desired normal abelian self-centralizing 
subgroup of G. 

Assume now that B ^ Cg{A). Since B is an abelian subgroup of G, it must be in list f|T3|) . 
So \B\ = p 2 , and either B = Z p x Z p or B = Z p 2 (the last case occurs only if p — 2). If 
B = Z p x Z p , then Cg(^4)/-B is a p-subgroup of GL 2 (p), in particular, |Cg(^4)/-B| = p. If 
B = Z p 2, then Cg{A)/B is a p-subgroup of Aut(Z, p 2). Since ^(p 2 ) = pQo — 1), it follows that 
\Cg(A)/B\ = p. So in any case we have |Cg(A)| = p 3 . 

Now the arguments depend on p. 

• In the case p = 2, we have that Cg{A) is a nonabelian group of order 8. Thus Cg{A) 
is either dihedral group D 8 or the quaternion group Q 8 . If Cg(A) is dihedral, then it 
possesses the unique (and hence characteristic) subgroup H = Z4, so if is the desired 
normal self-centralizing subgroup of G. If G = Qg is quaternion then, as we describe in 
Section 16.3.31 trying to impose a Q$ symmetry group on the 3HDM potential will result 
in a potential symmetric under a continuous group. Thus, this situation cannot happen 
if we search for finite realizable groups G. Note that this feature is purely calculational 
and does not rely on the existence of a normal maximal abelian subgroup which we prove 
here. 

• In the case p = 3, we have that Cg{A) is a nonabelian group of order 27 and exponent 
3, i.e. for every g G G we have g 3 = 1. In this case we do not yet know whether 
B is normal in G, but it is definitely normal in its own normalizer B < Ng{B) < G. 
Moreover Cq{A) < Ng(B), since B is normal in Cg{A). We can then consider the 
factor group Ng(B)/B. B = Z 3 x Z 3 is a maximal abelian group in PSU(3); therefore 
it is self-centralizing in PSU(3) and, consequently, in G. Then, in particular, we have 
that Ng(B)/B is a subgroup of Aut(B) = GL 2 (3). Moreover, the analysis which will 
be exposed in detail in Section [7] allows us to state that N PSU ^(B)/B = SX 2 (3), so 
Nq(B)/B is a subgroup of SX 2 (3). We show in Section [7| that one cannot use elements of 
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order 3 from SL 2 (3) because the potential will then become invariant under a continuous 
symmetry group. Therefore, Nq(B)/B cannot have elements of order 3, which implies 
B is a Sylow 3-subgroup of Nq(B). The same statement holds for every group that 
lies "between" Ng(B) and B, in particular, to Cg(A). This contradicts to the fact that 
\Cg(A) : B\ = 3 and Co (A) < Nq(B). So this case is impossible. 

Summarizing the group-theoretic part of our derivation, we proved that any finite group G 
which can be realized as a Higgs-family symmetry group in 3HDM is solvable, and in addition 
it contains a normal self-centralizing abelian subgroup A belonging to the list ( fl3l) . Then, 
according to (J9j) the group G can be constructed as an extension of A by a subgroup of 
Aut(A). 

This marks the end of the group-theoretic part of our analysis. We now need to check all 
the five candidates for A, whose explicit realization were already given in [21], and by means 
of direct calculations see which extension can work in 3HDM. 

5 Detecting continuous symmetries 

Before we embark on analyzing each particular abelian group and its extensions, let us discuss 
an important issue. In this paper, we focus on discrete symmetries of the scalar sector in 
3HDM. The symmetry groups we study must be realizable, that is, we need to prove that a 
potential symmetric under a finite group G is not symmetric under any larger group containing 
G. In particular, we must prove that a given G-symmetric potential does not have any 
continuous symmetry. 

In principle, it would be desirable to derive a basis-invariant criterion for existence or 
absence of a continuous symmetry. Such condition is known for 2HDM, [T91I20] . while for 
the more than two doublets a necessary and sufficient condition is still missing. However, in 
certain special but important cases it is possible to derive a sufficient condition for absence 
of any continuous symmetry. Since this method relies on the properties of the orbit space in 
3HDM, we start by briefly describing it. 

5.1 Orbit space in 3HDM 

The formalism of representing the space of electroweak-gauge orbits of Higgs fields via bilinears 
was first developed for 2HDM, [HlEOlEB] , and then generalized to iV doublets in [13]. Below 
we focus on the 3HDM case. 

The Higgs potential depends on the Higgs doublets via their gauge-invariant bilinear com- 
binations (p^b, a,b = 1, 2, 3. These bilinears can be organized into the following real scalar r$ 
and real vector n, i — 1, ... ,8: 

i^M + (40g) + (4<pa (<p\<Pi) - (4fa) (0I01) + (402) - 2(4^) 
ro ~ 7! ' rs " 2 ' rs " 27! ' 

ri = Re(0|0 2 ) , r 2 = Im(0l0 2 ) , r 4 = Re(0^i) , 

r 5 = Im(<40i) , r 6 = Re(0^ 3 ) , r 7 = lm(0^ 3 ) . (14) 
The last six components can be grouped into three "complex coordinates" : 

rn = (0i02) = n + ir 2 , r 45 = (0^0i) =r 4 + ir 5 , r 67 = (0J,0 3 ) = r 6 + ir 7 . (15) 
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It is also convenient to define the normalized coordinates = rj/r . The orbit space of the 
3HDM is then represented by an algebraic manifold lying in the 1 + 8-dimensional euclidean 
space of r and r 4 and is defined by the following (in)equalities, [13]: 

r > , n 2 < 1 , \Z3dijkninjrik = , (16) 

where d^ k is the fully symmetric SU(3) tensor. It can also be derived that |n| is bounded 
from below: 

n 2 = a, -<a<l. (17) 

The value of a parametrizes S"?7(3)-orbits inside the orbit space. In particular, we will use 
this relation below when substituting r 2 + r| by ar\ — |r 12 | 2 — | r- 45 1 2 — |r 67 | 2 . 

Any U(3) transformation in the space of doublets (pi, <p2, 4>3 leaves tq invariant and induces 
an SO(8) rotation of the vector r^. Note that this map is not surjective, namely not every 
50(8) rotation of can be induced by a U(3) transformation in the space of doublets. There- 
fore, unlike in 2HDM, we do not expect the orbit space of 3HDM to be S*0 (8)-symmetric, and 
the last condition in ( TT6"j) stresses that. 




P 



Figure 1: The orbit space of 3HDM in the (713, ng)-subspace (all other = 0). The outer and 
inner circles correspond to \n\ = 1 and \n\ = 1/2, respectively. 

Let us take a closer look at the (713, ns)-subspace. It follows from (ITBl that the orbit space 
intersects this plane along the equilateral triangle shown in Fig. CD Its vertices P, P', P" lie on 
the "neutral" manifold, which satisfy the condition n 2 = 1 and which would correspond to the 
neutral vacuum if the minimum of the potential were located there, while the line segments 
joining them correspond to the charge-breaking vacuum, see details in [13]. The orbit space in 
this plane clearly lacks the rotational symmetry and has only the symmetries of the equilateral 
triangle. 

5.2 Absence of continuous symmetries 

The convenience of the formalism of bilinears is that the most general Higgs potential becomes 
a quadratic form in this space: 

V = -M r - Min + ^Aooro + A 0l r r i + \Kf l r J . (18) 
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The real symmetric matrix Ay has eight real eigenvalues (counted with multiplicity). In order 
for the potential to be symmetric under a continuous group of transformations, Ay must have 
eigenvalues of multiplicities > 1. Note that any statement about eigenvalues of Ay is basis- 
invariant and therefore it can be checked in any basis. Furthermore, if we find a basis in 
which Ay has a block-diagonal form, and if eigenvalues from different blocks are distinct, then 
a continuous symmetry requires that each block is either invariant under this symmetry, or 
contains eigenvalues with multiplicity > 1. 

Let us consider an important special case of this situation. Suppose that the potential 
has no terms of type {<p\4>a){.4>b ( t ) c)^ where a,b,c are all distinct. This implies the absence of 
terms ?"o,3,8 r i, 2,4,5,6,7; an d the block- diagonal form of Ay, in which two blocks correspond to 
the (r3,r§) subspace and to its orthogonal complement. Suppose also that the eigenvalues of 
Ay in the (V^rg) subspace are distinct from those in the orthogonal complement. It follows 
then that any possible continuous symmetry must act trivially in the (r 3 , r 8 ) subspace, because 
the orbit space here lacks the rotational invariance. However, if ro, r%, and r§ are fixed, then 
4>i4>i, 02^2, and 0303 are also fixed. So, the doublets do not mix, and the possible continuous 
symmetry group can only be a subgroup of the group of pure phase rotations, which were 
studied in [21]. 

If in addition it is known that a given potential is not symmetric under continuous phase 
rotations, then we conclude that it does not have any continuous symmetry from PSU (3). 
It turns out that all the cases of various finite symmetry groups we consider below, except 
the last one, are of this type. Therefore, these arguments guarantee that the corresponding 
potentials can have only finite symmetry groups. Absence of a continuous symmetry in the 
very last case will be proved separately. 

6 Possible extensions: the torus chain 

We now check all the candidates for A from the list f[T3"j) and see which extension can work 
in 3HDM. In this section we will deal with the first four groups from the list, which arise as 
subgroups of the maximal torus; the last group will be considered later. For each group A, we 
use its explicit realization given in [21J as a group of rephasing transformations, and then we 
search for additional transformations from PSU (3) with the desired multiplication properties. 

6.1 Representing elements of PSU(S) 

Before we start analysis of each case, let us make a general remark on how we describe the 
elements of PSU (3). Using the bar notation for the canonical homomorphism SU(3) — > 
PSU '(3), we denote H < PSU (3) if its full preimage in SU(3) is H. Denoting the center of 
SU{3) as Z = Z{SU{3)) ~ Z 3 , we have Z = {1, z, z 2 }, where 

z = diag(w, co, to), u = e 2 ™ /3 . (19) 

The elements of the group H (a,b, . . . e H) will be written as 3 x 3 matrices from SU(3). The 
elements of H (a, b, ■ ■ ■ e H) are the corresponding cosets of Z in H . Explicit manipulation 
with these cosets is inconvenient, therefore in our calculation we represent an element a G 
PSU (3) by any of the three representing elements from SU(3): a, az, or az 2 . We will usually 
choose a and then prove that this representation is faithful (does not depend on the choice of 
representing element). 
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6.2 Extending Z 2 and Z 3 



The smallest group from the list is A = Z2, whose automorphism group is Aut{%2) = {1}, so 
that G = Z2. This case was already considered in [2T] . 

The next possibility is A = Z3, whose Aut(^L$) = Z2. The only non-trivial case to be 
considered is G/A = Z2, which implies that G can be either or — S3, the symmetry 
group of the equilateral triangle. The former can be disregarded because it does not appear 
in the list ({TBI , thus we focus only on the D 6 case. 

6.2.1 Constructing D 6 

The group D$ is generated by two elements a, b with the following relations: a 3 = 1, b 2 — 1, 
a& = 6a 2 . Following [21] . we represent the Z3 group by phase rotations: 

a = diag(w,w 2 ,l) . (20) 

There are in fact three such groups which differ only by the choice of the doublet which is fixed. 
However their generators, a, az, and az 2 , differ only by a transformation from the center, and 
therefore all of them correspond to the same generator a from PSU(3). It is straightforward 
to check that selecting a to represent a is a faithful representation. 

The explicit solution of the matrix equation ab = ba 2 shows that b G SU(3) must be of the 
form 

/ e iS \ 
6= e- iS , (21) 

V 0-1/ 

with arbitrary 8. The choice of the mixing pair of doublets (0i and 02 i n this case) is fixed by 
the choice of invariant doublet in a. 

The fact that b is not uniquely defined means that there exists not a single D e group but 
a whole family of Dq groups parametrized by the value of 8. Below, when checking whether 
a potential is D e symmetric, we will need to check its invariance under all possible D$s from 
this family. 

The generic Z3-symmetric potential contains the part invariant under any phase rotation 

with all indices running from 1 to 3, and the following additional terms 

Vi 3 = \i(4<h.)(4<l>i) + A 2 (0 t 1 02)(0^2) + + h.c. (22) 

with complex Ai, A 2 , A3. At least two of them must be non-zero, otherwise the potential will 
be symmetric under a continuous group of Higgs-family transformations, [2T] . Let us denote 
their phases by ipi, ip2, and ^3, respectively. If the parameters of Vo satisfy 

m ii = m 22 ) ^11 = ^22 , A13 = A23 , A' 13 = A 23 , (23) 

and if, in addition, |Ai| = IA2I, then the whole potential becomes symmetric under one partic- 
ular Dq group constructed with b in (|2ip with the value of 5 = (ip2 — tpi + 7r )/3 + 2iik/3. The 
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extra freedom given by 2irk/3 corresponds to three order-two elements of D e : b, ab, a 2 b. We 
opt to define b by setting k = 0. Alternatively, we can be compactly write the condition as 



35 = 7T - ipi + ip 2 ■ 



(24) 



To summarize, the criterion of the D e symmetry of the potential is that, after a possible 
doublet renumbering, conditions fl23|) and (124"]) are satisfied. 

Let us also note that when constructing the group D$ we could have searched for b satisfying 
not ab = ba 2 but ab = ba 2 ■ z p , with p = 1,2. Solutions of this equations exist, but they do 
not lead to any new possibilities. Indeed, let us introduce a' = az p . Then, we get a'b = ba' 2 . 
Thus, we get the same equation for b as before, up to a cyclic permutation of doublets, the 
possibility which we already took into account. 

6.2.2 Proving that D G is realizable 

This construction allows us to write down an example of the Dg-symmetric potential: it is Vq 
restricted by conditions ( 1231) with arbitrary choice of the doublet pair plus Vi 3 in ( 1221) subject 
to | Ai| = | As | . In order to show that _D 6 is realizable, we need to demonstrate that (at least 
in some cases) this potential is not symmetric under any larger Higgs-family transformation 
group. 

This proof is short and contains two steps. First, we note that the conditions described 
in section are fulfilled: the (r 3 , r 8 )-subspace does not couple to its orthogonal complement 
via Aij, and that the eigenvalues in these two subspaces are defined by different sets of free 
parameters. The extra terms (122]) guarantee that there is only finite group of phase rotations, 
the group Z 3 . Therefore, the generic Z) 6 -symmetric potential has no continuous symmetry. 

Second, we need to show that the generic Dg-symmetric potential has no higher discrete 
symmetries. This is proved by the simple observation that all other finite groups to be dis- 
cussed below which could possibly contain D§ lead to stronger restrictions on the potential 
than fl23|) and |Ai| = | As | - Therefore, not satisfying those stronger restrictions will yield a 
potential symmetric only under Dq. 

6.2.3 Including antiunitary transformations 

Any generalized-CP (antiunitary) transformation acting on three doublets is of the form 



Here J is the operation of hermitean conjugation of the doublets. If G is the symmetry group 
of unitary transformations, then it is normal in (G, J'), and J' induces automorphisms in G. 
So, when we search for J', we require 



where a generically denotes the generators of G. If such a transformation is found, the group 
is extended from G to G X 7T 2 , where asterisk on the group indicates that its generator is 
antiunitary. Note the crucial point of our method: when extending G by an antiunitary 
transformation, we require that the unitary transformation symmetry group remains G. 

Now, turning to extension of D 6 by an antiunitary symmetry, we first note that the resulting 
group Dq x is a non-abelian group of order 12 containing a normal subgroup D§. Among 



J' — c - J , ce PSU(3) . 



(25) 



{J'feG, (J')-WeG, 



(26) 



18 



the three non-abelian groups of order 12, there exists only one group, namely D 6 x Z^j, with 
a sugroup D G (which is automatically normal because all subgroups of index 2 are normal). 
This fact can also be proved in a more general way without knowing the list of groups of order 
12. Note that it contains, among other, the subgroup Zg; its presence does not contradict the 
list ( TT3"j) because that list refers only to the groups of unitary transformations. 

Next, let us denote the generator of Z^ by J' = cJ. Since J' centralizes the entire D 6 , it 
follows that (J')~ l aJ' = a, (J')~ 1 bJ' = b, and (J') 2 = cJcJ = cc* = 1. The matrix c satisfying 
these conditions must be of the form 

/ \ 

c = I , (27) 

\ -e~ 2i ~< / 

with arbitrary 7. Requiring the potential stays invariant under J', we obtain the following 
conditions on 7: 67 = + ^2) = 2ip 3 . Therefore, if the following extra condition is 

fulfilled, 

2(^1 + ^2 + ^3) = 0. (28) 

the .De-invariant potential becomes symmetric under the groups Dq x Z^. If this condition is not 
satisfied, the symmetry group remains Dq even in the case when antiunitary transformations 
are allowed. We conclude that both Dq and Dq x Z2 are realizable in 3HDM. 

It is interesting to note that if we set A3 = 0, then the potential would still be invariant 
under D§. However in this case it becomes symmetric under J' with 67 = —2(tpi + ^2), 
without any extra condition on ip\ and ip2, and the potential becomes automatically invariant 
under D§ x Zg. So, we conclude that the fact that Dq is still realizable even if anti-unitary 
transformations are included is due to the special role of the Z3-symmetry: we have three, not 
two terms in the Z 3 -symmetric potential, and it is the third term that prevents an automatic 
anti-unitary symmetry. 

6.3 Extending Z4 

Let us now take A = Z 4 (generated by a). Then Aut{7L$) = Z 2 , so that G = Z 4 . Z 2 generated 
by a and some b ^ Z 4 . The two non-abelian possibilities for G are the dihedral group D 8 
representing symmetries of the square, and the quaternion group Q%. In both cases b^ab = a 3 , 
with the only difference that b 2 = 1 for Dg while b 2 = a 2 for Q$. Note that extension leading 
to the dihedral group is split, D$ = Z 4 x Z 2 , while Q 8 is not. 

6.3.1 Constructing D 8 

Representing a by phase rotations a = diag(z, —i, 1), we find that b satisfying these conditions 
is again of the form (12 ip with arbitrary 5. However now we do not have the freedom to choose 
the pair of doublets which are mixed by b: this pair is fixed by a. Also, unlike the Z3 case, 
the matrix equation ab = ba 3 ■ z does not have solutions for b G SU(3). 
The Z 4 -symmetric potential (for this choice of a) is Vq + V^ 4 , where 

V %A = Xi(44>i)(4h) + H4h) 2 + h.c. (29) 
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The phases of Ai and A 2 are. as usual, denoted as ipx and ip 2 , respectively. Upon 6, the first 
term here remains invariant, while the second term transforms as 

(0l0 2 ) 2 ^ e-^(4<Pi) 2 - (30) 

This means that the potential fl29|) is always symmetric under (T2T]) provided that we choose 

5 = ^/2, (31) 

Therefore, in order to get a D 8 -symmetric potential we only require that Vo satisfies conditions 
(123]) . The proof that Dg is realizable (as long as only unitary transformations are concerned) 
follows along the same lines as in section 16.2.21 



6.3.2 Including antiunitary transformations 

In [21] we found that exactly the same conditions, namely ( 1231) and (T3"Tj) . must be satisfied 
for existence of an antiunitary transformation commuting with the elements of Z4. This 
transformation is again J' = cJ with c given by ( 12"T|) with 67 = 2ipi, and it commutes with 
all elements of Dg. Therefore, if we include antiunitary transformations, we automatically get 
the group Dg x Z^, while Dg becomes non-realizable. Note that the resulting group does not 
contain Zg. Indeed, we showed in [21] that imposing Zg symmetry group leads to a potential 
with continuous symmetry. 

6.3.3 Attempting at Q 8 

Solving matrix equations ah = ba 3 and b 2 = a 2 , we get the following form of b: 






e iS 


e -i5 








1 



b(Qg) = -e- lb 0. (32) 



By checking how V^ 4 in ([29]) transforms under it, we find that the first term simply changes 
its sign. The only way to make the potential symmetric under Qg is to set Ai = 0. But then 
we know from [21] that the potential becomes invariant under a continuous group of phase 
rotations. Therefore, Qg is not realizable. 



6.4 Extending Z 2 x Z 2 

If A = Z 2 x Z 2 , then Aut(Z 2 x Z 2 ) = GL 2 (2) = S3. The group Z 2 x Z 2 can be realized as 
the group of independent sign flips of the three doublets with generators a x = diag(l, —1, — 1) 
(equivalent to the sign flip of the first doublet) and a 2 = diag(— 1, 1, — 1) (equivalent to the 
sign flip of the second doublet), so that a\a 2 is equivalent to the sign flip of the third doublet. 
The potential symmetric under this group contains Vq and additional terms 

Vz 2 xz 2 = A 12 (0 t 1 2 ) 2 + A 23 (0^ 3 ) 2 + A 3 i(0 3 0i) 2 + h.c. (33) 

with at least two among coefficients Ajj being non-zero. The coefficients can be complex; as 
usual we denote their phases as ipij. This model is also known as the Weinberg's 3HDM, [5]. 
The non-abelian finite group G can be constructed as extension of A by Z 2 , by Z 3 or by £3. 
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6.4.1 Extention (Z 2 x Z 2 ) . Z 2 



Consider first the extension (Z 2 x Z 2 ) . Z 2 . The only extension leading to a non-abelian group 
is (Z 2 x Z 2 ) . Z 2 = D 8 , and we already proved that this group is realizable. Nevertheless, we 
prefer to explicitly work it out to see the reduction of free parameters. 

The element b which we search for must act on {ai,a 2 ,aia 2 } as a transposition of any 
pair. In addition, b 2 G Z 2 x Z 2 . It does not matter which pair of generators is transposed, as 
this choice can be changes by renumbering the doublets. So, we take b such that b _1 aib = a 2 
and b~ x a 2 b = a%. Then, b 2 can be either 1 or a\a 2 , because choices b 2 = a\ or a 2 lead to 
inconsistent relations. Indeed, if we assume b 2 — a\, then 

a 2 = b~ a\b = b~ x b 2 b = b 2 = a\ , 

which is a contradiction. In both cases (b 2 = 1 and b 2 = aia 2 ) we get the group D 8 . Even 
more, we get the same D 8 group: if b 2 = a\a 2 , then b' = bai satisfies b' 2 = 1, while its action 
on ai and a 2 remains the same. So, it is sufficient to focus on the b 2 = 1 case only. 

Again, explicitly solving the matrix equations, we get b of the form ( I2T1) with arbitrary 5. 
Then, we check how the potential (|33|) changes upon b and find that we need to set 

45 = 2Vi 2 , 26 = -(^23 + ^31), |A 23 | = |A 3 i|. (34) 

Equations on the phase 5 can be satisfied if 

2(^i2 + ^2 3 + ^3i) = & Im(\ 12 \ 23 \ 31 ) = 0. (35) 

So, if: (1) this condition is satisfied, (2) two among |Ay| are equal, (3) condition on Vo ( |23"j) is 
satisfied, then the potential is Dg-symmetric. Note also that if A12 = (which we are allowed 
to consider because ( 1331) contains three rather than two terms), then condition on the phases 
is not needed. 

It might seem that these conditions on the potential to make it Dg-symmetric are more 
restrictive than in the Z 4 extention we studied above. However note that the Z 2 x Z 2 -symmetric 
potential (|33|) has six free parameteres, and we placed two conditions to reduce the number 
of free parameters in the Dg potential to four (apart from Vo). On the other hand, ( )29l) had 
only four from the beginning, and without any restriction this number survives. Therefore we 
have the same number of degrees of freedom when constructing D 8 in either way. 



6.4.2 Constructing (Z 2 x Z 2 ) x Z 3 = T 

The extension by Z 3 is necessarily split, (Z 2 x Z 2 ) x Z 3 , leading to the group T ~ A±, the 
symmetry group of the tetrahedron. To construct it, we need b such that b 3 = 1 with the 
property that b acts on {ai, a 2 , a\a 2 } by cyclic permutations. Fixing the order of permutations 
by b^aib = a 2 , we find that b must be of the form 

/ e iSl \ 
b = I e l5 * , (36) 

\ e ^ Sl+5 ^ / 

with arbitrary 5%, 5 2 . It then follows that if coefficients in f )33|) satisfy 

|Ai 2 | = I A23 1 = I A31 1 , (37) 
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then Vz 2 xz 2 is symmetric under one particular b with 



2^12 - ^31 - ^23 r 2^23 - V>31 ~ ^12 

<h = ^ , o 2 = 



6 



6 



Then, by a rephasing transformation one also make the phases of all \j equal and bring (1331) 
to the following form 



V T = X 



>I0 2 ) 2 + (4fo) 2 + (4<h) 2 



+ h.c. 



(38) 



with a complex A. In this form, the parameters S± = 82 — 0, and the matrix b is just the cyclic 
permutation of the doublets. In addition, the symmetry under b places stronger conditions on 
the parameters of Vo, so that the most general Vq satisfying them is 



Vo 



-m 



10iJ + 



402) 



+ 



403) 



+ A 



i0ij + 



402) 



+ 



403) 



(39) 



+A' 



(0l0l)(0|0 2 ) + (0^02)(0|03) + (0303) (0l0l)] + A" (|<M0 2 | 2 + |0 2 3 | 2 + |0 3 



6.4.3 Constructing (Z 2 x Z 2 ) x S 3 = O 

The last extension, (Z 2 x Z 2 ) . ^3, is also split, otherwise we would obtain Zg. It leads to the 
group O ~ S4, the symmetry group of the octahedron and the cube. As it includes T as a 
subgroup, the most general O-symmetric potential is Vq from (!39l) plus Vr from (J38l) with the 
additional condition that A is real (the extra symmetry with respect to the T-symmetric case 
is a transposition of any two doublets). 



6.4.4 Including antiunitary transformations 

The case of D$ has been already considered in section I6.3.2L 

The tetrahedral potential Vr + V from (1381) and fl39|) is symmetric under the following 
antiunitary transformation: 

/0 1 \ 

J' = 1 J, (40) 

\o -1/ 

which generates a Zg group. Therefore the symmetry group of this potential is the full achiral 
tetrahedral group ~ T x Zg, which is isomorphic to S4. 

The octahedral potential is a particular case of the tetrahedral one, therefore it is also 
invariant under an antiunitary transformation. The extra Z 2 subgroup is generated by the 
complex conjugation, J, and this transformation commutes with the entire Higgs-family group 
O. Therefore, the symmetry group of the potential is the full achiral octahedral symmetry 
group O h ~ O x Z^. 



6.5 Extensions of abelian groups by an antiunitary transformation 

The last type of extension we need to consider is of the type A . Z 2 , where A is one of the four 
abelian groups of Higgs-family transformations lying in a maximal torus, that is, the first four 
groups in the list f)13p . while the Z?; is as usual generated by an antiunitary transformation 
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J' = cJ. This problem was partly solved in |21J, where such extensions leading to abelian 
groups were analyzed. It was established that only the following four abelian groups of this 
type are realizable: Zj, Z|, Z 2 x Zg, and Z 2 x Z 2 x Z^. Here, we consider non-abelian extensins 
of this type. 

6.5.1 Anti-unitary extension of Z 3 

The smallest non-abelian group we can have is Z3 x Z 2 — Dq. We stress that this D§ group 
we search for is different from what we analyzed in section I6.2[ because there the D e group 
contained only unitary transformations. Using the same notation for the generator a of the 
Z 3 group, we find that the transformation c in the definition of J' must be diagonal: c = 
diag(e 4 ^ 1 , e^ 2 , e -l ^ 1+ ^). Then, studying how the Z3-symmetric potential Vq + Vz 3 changes 
under J' = cJ, we obtain that the only condition to be satisfied is (12 8j) . 

If this condition is satisfied, then the potential is invariant under Z3 x Z 2 ~ Dq, if not, 
then the symmetry group remains Z 3 . This proves that both groups are realizable in 3HDM. 
Note that in contrast with the Dq x Z2 case, we do not place any extra condition such as (123]) . 

6.5.2 Anti-unitary extension of Z 4 

A priori, the two non-abelian extensions here are again Dg and Qs- With the usual convention 
for a, the generator of Z4, we again obtain that c must be of the same diagonal form. This 
immediately excludes the Qs case because we have (J') 2 = c*c = 1. 

The case of Z 4 x Z 2 ~ D 8 is possible. Even more, it turns out that the Z/t-symmetric 
potential V + Vz 4 is always symmetric under some J' of this type. It means, therefore, that 
if anti-unitary transformations are included, Z4 is not realizable anymore: the true symmetry 
group of the potential is Z4 x Z 2 ~ D$. In more physical terms, we conclude that presence of 
a Z4 group of Higgs-family transformations makes the potential explicitly CP -conserving. 

6.5.3 Anti-unitary extension of Z 2 x Z 2 

The only non-abelian extension of the type (Z 2 x Z 2 ) . can produce only D 8 , which was 
already considered. We only remark here that c turns out to be of the type (T2"T|) . which places 
extra contraints on Vq. Not satisfying these constraints will keep the symmetry group Z 2 x Z 2 , 
which means that it is realizable. 

7 The Z 3 x Z 3 chain 

7.1 The group and its extensions 

The last abelian group from the list ( !T3l) . Z3 x Z3, requires a special treatment due to a number 
of reasons. First, it does not belong to any maximal torus of PSU(3) but is a maximal abelian 
subgroup of PSU(3) on its own, [21], and its full preimage in SU(3) is the non-abelian group 
A(27), [29]. Second, its automorphism group Aut(Z 3 x Z 3 ) is sufficiently large and requires 
an accurate description. 
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Let us first remind how this group is constructed. We first consider the subgroup of SU (3) 
generated by 



: 
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1 


: 









!■ H 










V o 





- 2 ) 




V 1 





o / 



This group known as A(27) is non-abelian because a and b do not commute, but their com- 
mutator lies in the center of SU(3): 

[a, b] = aba^b- 1 = z 2 G Z{SU(3)) . (42) 

Therefore, its image under the canonical homomorphism SU(3) — > PSU(3) becomes the 
desired abelian group A(27)/Z 3 = Z 3 x Z 3 . The true generators of Z 3 x Z 3 are cosets a = 
aZ(SU{3)) and b = bZ(SU{3)) from PSU(3), and they obviously commute: [a, b] = 1. Note 
that since Z 3 x Z 3 is a maximal abelian subgroup in PSU(3), there is no other element in 
PSU(3) commuting with all elements of this group, so Cpsu(3)(^3 x Z 3 ) = Z 3 x Z 3 . 

If the normal self-centralizing abelian subgroup of G, whose existence was proved in sec- 
tion I4.2[ is A = Z 3 x Z 3 , then G can be constructed as an extension of A by a subgroup of 
Aut(Z 3 x Z 3 ) = GL2(3), the general linear group of transformations of two-dimensional vector 
space over the finite field F 3 . The order of this group is |GL 2 (3)| = 48, and it will prove useful 
if we now digress and describe the structure of this group in some detail. 



7.1.1 Z 3 x Z 3 as a vector space over F 3 

The finite field F 3 is defined as the additive group of integers mod 3, in which the multiplication 
is also introduced. It is convenient to denote the elements of this field as 0, 1, —1 with obvious 
addition and multiplication laws. Unlike the integers themselves, F 3 is closed under division 
by a non-zero number, the property that makes F 3 a field. 

A vector space over a finite field is defined just as over any "usual" field. The group Z 3 x Z 3 
can be thought of as a 2D vector space over F 3 ; its elements are (with the additive notation 
for the group operation) x = q a a + q b b, where q a , q b G F 3 , and a, b are, as before, the generators 
of the group Z 3 x Z 3 . In the multiplicative notation, we write x = a qa b qb . 

It is possible to define an antisymmetric scalar product in this space. For any x G Z 3 x Z 3 , 
take any element of its preimage, x G A(27). Then, for any two elements x, y G Z 3 x Z 3 , 
construct the number (x,y) as [x, y] G F 3 . This map is faithful: although we can select 
different x for a given x, all of them give the same [x, y]. 

Clearly, (x,y) = —(y,x), in the additive notation. Besides, the so defined product is linear 
in both arguments: 

(x x + x 2 ,y) = (x u y) + (x 2 ,y) , (x, y x + y 2 ) = (x, yx) + (x, y 2 ) ■ (43) 
Indeed, for any three elements of any group the following relation holds: 

[xy, z] = xyzy~ 1 x~ 1 z~ 1 = xyzy^ 1 ■ z~ l x~ l xz • x^ 1 z^ 1 = x[y, z]x~ x [x, z] . (44) 

If in addition all commutators take values in the center of the group SU(N), then x and x _1 
can be cancelled, and we get [xy, z] = [y, z\[x, z\. In our case we represent x x = a gal b qbl z ri and 
similarly for x 2 and y, and noting that all z ri are inessential, we recover the above linearity 
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in the first argument. Thus, Z 3 x Z 3 becomes a vector space over F 3 equipped with an 
antisymmetric scalar product. 

Note that all antisymmetric products in Z 3 x Z 3 are proportional to (a, b). Indeed, if two 
elements x and x' are defined by their vectors q = (q a , g&) and <f = (q' a , q' b ), then due to 
bilinearity we get 

(x, x') = (q a q' b - q b q' a )(a, b) = e^q^a, b) , (45) 
where is the standard antisymmetric tensor with e\ 2 = — e 21 — 1, e±i — e 22 = 0. 



7.1.2 The automorphism group of Z 3 x Z 3 

The automorphism group of Z 3 x Z 3 can then be viewed as the group of non-degenerate 
matrices with elements from F 3 acting in this 2D space, which explains why Aut(Z 3 x Z 3 ) = 
GL 2 {3). Each matrix q can be defined by its action on the generators a, b: o4 q aa a + q a bb, 
b i — y qbofl + qbbb, and can therefore be written as 

q=( qaa qab ) , detg^O. (46) 

\ Qba qbb ) 

The group operation in GL 2 {3) is just the matrix product. 

Recall now that the elements of both the Z 3 x Z 3 group and of its automorphism group are 
represented in our case as unitary or antiunitary transformations of the three doublets (that is, 
we work not with the abstract groups but with their 3-dimensional complex representations). 
Since Z 3 x Z 3 is assumed to be normal in G, the elements g G Aut(Z 3 x Z 3 ) act on the elements 
of Z 3 x Z 3 by conjugation: x h- > g~ x xg, which we denoted by g{x). Then the antisymmetric 
product defined above changes upon this action in the following way: 

(g(xU(V)) = g- 1 ^ v]g = g~^ r g = { [/y t^y = _ (i> v) , !^ I ZtZLry . 

(47) 

Here we used the fact the commutator of any two elements of A(27) lies in the center Z(SU (3)), 
and that the CP conjugation operator J acts on any x G SU(3) by J~ x xJ = x*. So, unitary 
transformations preserve the antisymmetric product, while anti-unitary ones flip its sign. 

Generically, the subgroup of a general linear group which conserves an antisymmetric 
bilinear product in a vector space is called symplectic. Here we have the group Sp2{3) < 
GL 2 (3). It turns out that Sp2{3) = SL 2 (3). Indeed, suppose g G GL 2 (3) acts in the 2D space 
over F 3 by mapping q^ i— )■ g(q) = gwqi'. Then, the product transforms as 

(x, y) H> (g(x),g(y)) = e i:j g iV g jr qPq^\a, b) = det g ■ (x, y) . (48) 

Since detg = ±1, we get two kinds of transformations: those which conserve all products 
(detg = 1, so that g G SX 2 (3)) and those which flip their signs (det g = —1), hence the 
identification of Sp2(3) and SX 2 (3) follows. 

We conclude that the finite symmetry group G with the normal self-centralizing abelian 
subgroup Z 3 x Z 3 can be constructed by extension (Z 3 x Z 3 ) . K, where K < SL 2 (3). 
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7.1.3 Explicit description of SL 2 (3) 



The structure of the group SL,2(3) is well-known, but it will prove useful later to have the 
explicit expressions for some of its elements. 

The order of the group is |SX 2 (3)| = 24. It contains elements of order 2, 3, 4, and 6, 
generating the corresponding cyclic subgroups. The subgroup Z 2 is generated by the center 
of the group 

-1 
-1 



(49) 



which in the multiplicative notation means a i— >■ a 2 , b y b 2 . There are four distinct Z 3 
subgrops generated by 



h ~ [I 1 ) ' h ~ 
three Z 4 subgroups generated by 

ll< ~ ^ 1 ^ 



1 

1 1 



h 



1 1 
1 -1 



1 
-1 -1 



-1 1 
1 1 



-1 

1 -1 



(50) 



(51) 



and four subgroups, which we do not write explicitly because they are absent in the list 

Each such element can be represented by a unique (up to center) SU(3) matrix, which can 
be found by explicitly solving the corresponding matrix equations defining the action of this 
element. For example, the transformation c is defined by 



z(a) = c 1 ac = a 2 , c{b) = c 1 bc = b 2 



(52) 



Rewriting these equations as 3 x 3 matrix equations ac = ca 2 , be = cb 2 and solving them 
explicitly, we find the matrix c: 



(53) 




7.2 Generic potential 

A generic potential symmetric under Z 3 x Z 3 is 



V 



—m 

+4 

+A 2 
+A 3 



>l0l) + (0202) + (0303) 



+ A c 



^0l) + (0^02) + (0^03) 



1 2 



(0l0l) 2 + (0202) 2 + (0^03) 2 - (0 f 10l)(0202) - (0202X0303) " (4fo)(4<l>l) 



(0 t 1 2 )(0 t 1 3 ) + (0203)(020l) + (0300(0302 



+ h.c. 



(54) 



with real m 2 , Ao, Ai, A 2 and complex A 3 . All values here are generic. This potential can 
be found by taking the potential symmetric under the Z 3 group of phase rotations described 
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above and then requiring that it be invariant under the cyclic permutations on the doublets. 
Written in the space of bilinears, the potential has the form 

V = - v ^ m 2 ro + 3Aor 2 + v / 3A 1 (r 3 2 + r 8 2 ) + A 2 (|r 12 | 2 + |r45| 2 + |r 67 | 2 ) 
+A 3 (r 12 r* 5 + r 67 r* 2 + r 45 r* 7 ) + A3 (r* 2 r 45 + r* 7 r 12 + r* 5 r 67 ) 
= -\/3m 2 r + 3\ rl + AijTiTj . (55) 

It is important to prove that this potential has no continous symmetry. Using the approach 
described in section [5j we calculate the eigenvalues of Ay and find that it has four distinct 
eigenvalues of multiplicity two: 

v^Ai, A 2 + A 3 + A;, X 2 + lo\ 3 + lo 2 \* 3 , X 2 + lo 2 \ 3 + lo\* 3 . (56) 

The first eigenvalue corresponds to the subspace (r 3 , r 8 ), while the rest are three 2D subspaces 
within its orthogonal complement (?"i, t 2 , T4, T5, T6> r 7)- F° r generic values of the coefficients 
they do not coincide. Then, according to our discussion in section a continuous symmetry 
group, if present, must consist only of phase rotations of the doublets. But the A 3 term selects 
only the Z 3 group of phase rotations, which proves that no continuous symmetry leaves this 
potential invariant. 

7.3 Extension (Z 3 x Z 3 ) x Z 2 

It turns out that Z 3 x Z 3 is not realizable because the potential ( 154p is symmetric under a 
larger group (Z 3 x Z 3 ) x Z 2 = A(54)/Z 3 , which is generated by a,b,c with the following 
relations 

a 3 = b 3 = 1, c 2 = 1, [a, b] = 1, cdc = a 2 , cbc = b 2 . 

In terms of the explicit transformation laws, c is the coset cZ(SU(3)), with c being the 
exchange of any two doublets, for example ( |53|) . Note that (a, c) = S3 is the group of arbitrary 
permutations of the three doublets. Thus, if G = (Z 3 x Z 3 ) . K, then a G-symmetric potential 
must be a restriction of (1541) . and K must contain a Z 2 subgroup. 

There are three kinds of subgroups of SX 2 (3) containing Z 2 but not containing Z 6 : Z 2 , 
Z 4 , and Qs- In each case it would give a split extension, so G must contain a subgroup 
isomorphic to one of these groups. Since, as we argued above, the quaternion group Q 8 is not 
realizable in 3HDM, K can only be Z 2 or Z 4 . Therefore, the only additional case to consider 
is (Z 3 x Z 3 ) xi Z 4 , the group also known as S(36), [29] . 

7.4 Extension (Z 3 x Z 3 ) x Z 4 

There are three distinct Z 4 subgroups in SX 2 (3) generated by di, <i 2 , and d 3 , listed in fl5T| . In 
principle, all of them are conjugate inside SX 2 (3), but for our purposes all of them need to be 
checked. Explicit solutions of the matrix equations give the following transformations: 

i { 1 1 1 \ if 1 1 U \ z/ lla;2 

d ± = — \ 1 u 2 oj , d 2 — —j= \ 00 1 1 , d 3 = —= \ u 1 oj 

^ 3 \ 1 OJ OJ 2 J ^ 3 \ oj oj 2 oj J ^ 3 \ 1 OJ OJ 

Note that the prefactor z/v3 can also be written as 1/(oj 2 — oj). 
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Let us mention here that when searching for explicit SU(3) realizations of the transfor- 
mations di, we solve equations d^adi = b, d^bdi = a 2 . However, we could also use other 
representative matrices, a' and b', which differ from a and b by transformations from the center. 
For example, we can also ask for solutions of 

d'-Wi = z ni b , d'^bd^ = z n2 a 2 . (58) 

However, the solution of this equation can be written as 

d[ = d ia ni b n2 . (59) 

Therefore the resulting group {d[, a, b) coincides with (di, a, b). The similar results hold for di 
and o?3. 

7.4.1 Conditions for the (Z 3 x Z 3 ) X Z 4 symmetry 

We should now check how the potential ( 154]) changes under these transformations and when 
it remains invariant. The calculation is simplified if we introduce the folowing combinations 
of bilinears (here i*j stands for 4>\4>j): 

A = 1*1 + 2*2 + 3*3, Ai = 1*1 + w2*2 + w 2 3*3 , A 2 = A\ 

B = 1*2 + 2*3 + 3*1 , B 1 = l*2 + w2*3 + w 2 3*l, B 2 = 1*2 + w 2 2*3 + u3*l , 

B* = 2*1 + 3*2 + 1*3, B{ = 2*1 + w 2 3*2 + ul*3 , B* 2 = 2*1 + w3*2 + u 2 l*3 . 

Next, introducing 



X 



1 

1 l^l 2 



V3 

y = + \<f>w + i<+'i 



U j2 , uU_|2 , uU .12 _ l^o| 2 + \Bi\ 2 + \B 2 \ 2 



z* = (<p\fcMfo) + (4foM<Pi) + tfMitfM = W + Uj2lB f + ujW , (60) 



we write the potential ( 154"]) as 

V = -V3m 2 r + 3A r 2 + \*X { , where A*JQ = X t X + A 2 F + \ 3 Z* + A*Z (61) 

is the scalar product of the vector of coefficients and the vector of coordinates. Now, it follows 
from explicit calculations that the action of di can be compactly represented by the following 
transformations: 



d\ : 


A 1 


-+B , 


Bo 


-> A\ , 


B l w 2 5 2 , 


5 2 ->• B{ 


d 2 : 


A x 




B 1 


->■ ujA\ , 




B 2 — >• i?o 


d 3 : 


A 1 


^B 2 , 


B 2 


-)■ uA\ , 


-Bo Bi , 
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or even more compactly 



d 3 



: \A, 


2 


•H- 


|B 


2 

) 


1^1 


2 


•H- 


l^ 2 | 


: \A ± 


2 




\Bl 


2 

) 




2 






: |^i 


2 




\B 2 


2 

5 


\B 


2 


•H- 


l#l| 



(62) 



Therefore, their action in the space of (X, Y, Z, Z*) is given by the following hermitean and 
unitary matrices 



T(di 



/ o Vz Vz Vz\ 

V3 2 -1-1 

\/3 -1 -1 2 

V V3 -1 2 -1 / 



/ o Vs 



T{d 2 



C0 2 y/3 U)\/3\ 



V3 
coy/3 



\oo 2 VS -co 2 



2 -co 2 
-co -1 
2co 



-co 
2co 2 
-1 / 



and T(d 3 ) = [T(rf 2 )]*. It can be also noted that T(<i 2 ) acts in the space of (X, Y, co 2 Z, coZ*) 
by the matrix T(d\). So, T(di), T(d 2 ) and T(d 3 ) represent the same type of transformations 
acting in the spaces (X, Y, Z, Z*), (X, Y, co 2 Z, coZ*), or (X, Y, coZ, co 2 Z*), respectively That 
is, if (x, y, z, z*) is an eigenvector of T(di), then (x, y, coz, co 2 z*) is an eigenvector of T(d 2 ) and 
(x, y, co 2 z, coz*) is an eigenvector of T(d 3 ). This observation restores the expected symmetry 
among the three types of Z 4 subgroups inside SL 2 (3). 

Since these matrices are hermitian and unitary, they act by pure reflections. If we want 
the potential to be symmetric under one of these di, it must induce the same transformations 
in the space of A, = (Ai, X 2 , A3, A3). Therefore, in order to find conditions that the potential is 
invariant under di, we need to find eigenvectors of T(di) corresponding to the eigenvalue —1 
and require that Aj's projection on these eigenvectors is zero. 

Consider first T{d\). It has two eigenvectors corresponding to A = —1: (—y/3, 1, 1, 1) and 
(0,0, 1, —1). Therefore, we obtain the following condition for the potential to be symmetric 
under d\\ 

Aa . (63) 



Similarly, for d 2 we have 



A3 is real and As 



co\ 3 is real and oo\i 



\/3Ai — A2 



(64) 



For d 3 we have the complex conjugate condition. Together, they merge into a single condition 
for the potential (I54p to be symmetric under (Z3 x Z3) x Z4: 



I A3 



\V3\i- A 2 



and 3 arg A3 





7T 



if v^Ai > A 2 
if v^Ai < A 2 



(65) 



Let us also mention that when these conditions are taken into account, the spectrum of 
the matrix given in ( 15 6 p becomes even more degenerate: it contains two eigenvalues of 
multiplicity four (we refer to this spectrum as 4 + 4). 



7.4.2 Absence of a continuous symmetry 

In order for the group (Z 3 x Z 3 ) x Z 4 to be realizable, we need to show that the potential ()5 
with parameters satisfying (16"5]) is not symmetric under any continuous group. 
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We first note that even if such a continuous symmetry group existed, it could only be U(l). 
Indeed, the spectrum of in our case is 4 + 4, while for U(l) x {7(1) and SU(2) it must be 
6 + 2, and for 50(3) it must be 5 + 3. 

Let us now consider, for example, the <ii-symmetric potential. Using Y^=i r i = ar oi where 
1/4 < a < 1 parametrizes S{7(3)-orbits in the orbit space, we can rewrite it as 

V = - v / 3m 2 r +(3A +av / 3Ai)rg- ^— ^ — — (ki2 - r 45 \ 2 + |r 45 - r 67 | 2 + |r 67 - ri 2 | 2 ) . (66) 

Suppose the potential ( 1661) is invariant under a {7(1) group of transformations of doublets, 
generated by the generator t from the algebra su(3). Since the potential ( 16 6 p is invariant 
under the S 3 group of arbitrary permutations of the doublets, then the same potential must 
be also invariant under other {7(1) subgroups which are generated by various t 9 , which are 
obtained by acting on T by g G If t 7^ t 9 (or to be more accurate, if their corresponding 
U(l) groups are different), then the continuous symmetry group immediately becomes larger 
than i7(l), which is impossible. Therefore, t 9 must be equal (up to sign) to t for all g G S3. 
In other words, S3 must stabilize the {7(1) symmetry group. 

There exist only two elements in the algebra su(3) with this property: 



and t 2 = 1 1 . (67) 





ti generates pure phase rotations. It is explicitly S3-invariant, therefore the corresponding 
{7(1) group is also invariant. t\ induces SO (3) rotations of the doublets around the axis 
(1,1,1). It is Z3-invariant, while reflections from S3 flip the sign of t\. However the {7(1) 
group is still invariant. Since t\ and t 2 realize different representations of S3, one cannot take 
their linear combinations. So, the list of possibilities is restricted only to t\ and ti themselves. 
The eigenvectors and eigenvalues of t\ are 

A = : [ 1 ] , A = v / 3:(w 2 |, A = -V3 : [ u \ . (6£ 

CO J \co 2 

The presence of zero eigenvalue implies that the combination <f)\ + 4>2 + 4>3 is invariant under 
the {7(1) group generated by t\, and so are 




n 



+ 02 + 03| , 101+^02 + ^31, |01+W0 2 + ^0 3 | , (69) 



which simply means that r\ + r 4 + r 6 and r 2 + r 5 + r 7 are, separately, invariant. So, if the 
potential depends only on r and these two combinations, then it is symmetric under the {7(1) 
generated by t 1 . The point is that our potential (I6"6"j) cannot be written via these combinations 
only, therefore it is not invariant under this group. 
Consider now t 2 . Its eigensystem is 



A = 2 : 1 , A = -1 : 1 and -1 . (70) 
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There is no zero eigenvalue, therefore no linear combination of 0's is invariant. The indepen- 
dent quadratic combinations are 

|</>l + 02 + 3 | 2 , 102 -0 3 | 2 , I20X-02-03I 2 , (^-4)(20 1 -0 2 -0 3 ). (71) 

In addition, there exists a triple product of 0's which is also invariant but it is irrelevant for 
our analysis because our potential contains only two 0's and two 0''s. These invariants can 
also be rewritten as the following linearly independent invariants (here pi = (fi\4>i): 

pi + 2r 6 , p 2 + 2r 4 , p 3 + 2r x , r 2 + r 5 + r 7 . (72) 

Despite the fact that we now have more invariants than in the previous case, it is still impossible 
to express ( |66i) via these combinations. This means that ( |66i) is not symmetric under t 2 . 

This completes the proof that the potential ( 154"1) subject to conditions ( |65i) is not invariant 
under any continuous group. 



7.4.3 Absence of a larger finite symmetry group 

Although the group-theoretic arguments guarantee that no other extension can be used, it is 
still instructive to check what happens if we try to impose invariance under other subgroups 
of SL 2 (3). 

Let us first note that if we try to impose simultaneous invariance under two among di 
(trying to get Q 8 ), we must set A 3 = 0. But then the potential has an obvious continuous 
symmetry, and our attempt fails. 

Next, let us assume that the potential is invariant under (Z 3 x Z 3 ) x Z 3 , where the last Z 3 
is generated by one of the generators / in (15"U|) . for example f — fi- Its representative matrix 
in SU(3) is 

1 to 2 1 
1 1 uj 2 
< 2 1 1 



./ 




1 . 



(73) 



An analysis similar to what was described above allows us to find the corresponding transfor- 
mation matrix in the space of X, Y, Z, Z*\ 



T(fi 





( 


v/3 


V3uj 2 


VSuj \ 


1 


V3 


2 


-UJ 2 


— UJ 


3 


y/3uj 2 


-u 2 


—UJ 


2 




V V3u 


—u 


2 


-uj 2 J 



(74) 



It leads to the following conditions for the potential to be symmetric under (Z 3 x Z 3 ) x Z 3 : 

A 2 — A 3 



A, 



An and Ai 



In the space of bilinears, the potential can then be compactly written as 



V = -V3m 2 r + (3A + v^Aia)^ + A 3 |r i2 + r 45 + 



r&7 



(75) 



(76) 
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The spectrum of Ay becomes of the type 6 + 2. This high symmetry hints at existence of a 
possible continuous symmetry of the potential, and it is indeed the case. For example, the 
following SO (2) rotations among three doublets, 4> a i— > R a b(a)(p{ ) , leave r 12 + r 45 + r 67 invariant: 

j / l + 2cosa l + 2cosa" l + 2cosa' \ 
R(a) = -l l + 2cosa' l + 2cosa l + 2cosa" , (77) 
\ 1 + 2 cos a" 1 + 2 cos a' 1 + 2 cos a ) 

with a G [0, 27r) and a' = a + 27r/3 and a" = a + 47r/3. Note that at a = 0, 27r/3 and Att/3 
we recover the Z 3 group (b). 

We conclude therefore that imposing invariance under Z 3 < SX 2 (3) makes the potential 
symmetric under a continuous group. In this way, we completely exhausted possibilities offered 
by SL 2 (3). 

7.5 Anti-unitary transformations 

We showed in section [7TT1 that antiunitary transformations correspond to elements of GL 2 (3) 
not lying in SL 2 (3) as they have negative determinant and flip the sign of the antisymmetric 
scalar product in A = Z 3 x Z 3 . The complex conjugation operator, J, acts in A by sending a 
to a 2 and leaving b invariant. Therefore, the corresponding matrix is 

Since any antiunitary transformation can be written as J' = qJ, where q is unitary, it follows 
that q must belong to SX 2 (3). 

Next, we need to find which q's can be used. Clearly, (J') 2 = qJqJ = qq* G SX 2 (3). If we 
are looking for an antiunitary symmetry of a (Z3 x Z3) x Z2-symmetric potential, then qq* 
must be either 1 or c, which generates the center of SL 2 (3). 

Let us first consider the second possibility. 

■«=(!!)■ then «*=(-, ?)■ < 79 > 

Using this to solve gg* = c, we get six possible solutions, but all of them have det q = — 1, 
that is, they do not belong to SL 2 (3). Therefore, the only possibility is qq* = 1. 

But then we can apply the results of our search for antiunitary transformations for the 
D e case. Our group (Z 3 x Z 3 ) x Z 2 contains the D 6 subgroup with 5 = ir. Therefore, we 
arrive at the conclusion: in order for our potential to be symmetric under an antiunitary 
transformation, we must require 

6 arg A 3 = . (80) 

If this criterion is satisfied, the symmetry group becomes (Z 3 x Z 3 ) x (Z 2 x Z 2 ); otherwise the 
group remains (Z 3 x Z 3 ) x Z 2 . Therefore, both groups are realizable in 3HDM. 

Now, consider the case of the extended symmetry group, (Z 3 x Z 3 ) x Z4 ~ £(36). In 
this case ( |80l) is satisfied automatically due to ( 1631) . We then conclude that in this case the 
realizable symmetry is £(36) x Z 2 . 
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8 Summary and discussion 



8.1 List of realizable finite symmetry groups in 3HDM 

Bringing together the results of the search for the abelian symmetry groups, [21], and of the 
present work, we can finally give the list of finite groups which can appear as the symmetry 
groups of the scalar sector in 3HDM. If only Higgs-family transformations are concerned, the 
realizable finite groups are 

Z 2 , Z 3 , Z 4 , Z 2 xZ 2 , As, D s , T~A 4 , 0~S 4 , 

(Z 3 x Z 3 ) x Z 2 ~ A(54)/Z 3 , (Z 3 x Z 3 ) x Z 4 ~ £(36) . (81) 

This list is complete: trying to impose any other finite symmetry group of Higgs-family trans- 
formations leads to the potential with a continuous symmetry. 

If both unitary (Higgs-family) and antiunitary (generalized-CP) transformations are al- 
lowed, the list becomes 

Z 2 , Z 3 , Z 2 x Z 2 , Z 2 , Z 4 , 

Z 2 x Z* 2 , Z 2 x Z 2 x Z\, Z 3 x Z* 2 ~ P> 6 , Z 4 x Z^ ~ P> 8 , 

D e , D 6 x Z* 2 , D 8 x Z* 2 , A 4 xZ* 2 ~T d , S±xZ*~O h , 

(Z 3 x Z 3 ) x Z 2 , (Z 3 x Z 3 ) x (Z 2 x Z 2 ), £(36) x Z* . (82) 

As usual, an asterisk here indicates that the generator of the corresponding group is an anti- 
unitary transformation. Note that Higgs-family transformation groups Z 4 , D 8 , A±, S 4 , and 
£(36) become non- realizable in this case, because potentials symmetric under them are au- 
tomatically symmetric under an additional anti-unitary transformation. In all cases apart 
from A A this is the consequence of our finding in section 16.51 that presence of the Z 4 group of 
Higgs-family transformations always leads to an additional anti-unitary symmetry. 

These lists complete the classification of realizable finite symmetry groups of the scalar 
sector of 3HDM. Conditions for the existence and examples of the potentials symmetric under 
each of these groups have been given in [21] and in the present work. For the reader's con- 
venience, we collected these examples involving non- abelian Higgs-family symmetry groups in 
the Appendix. 

8.2 Interplay between Higgs-family symmetries and explicit CP- 
violation 

In 2HDM, presence of any Higgs-family symmetry immediately implies presence of a generalized- 
CP symmetry. In other words, it is impossible to write down an explicitly CP-violating 2HDM 
potential with any Higgs-family symmetry. In this sense, generalized-CP symmetries can be 
viewed as the smallest building blocks of any symmetry group in 2HDM. 

By comparing lists ( 18T1) and ( |82|) . we see that this conclusion is no longer true in 3HDM, 
namely there are some Higgs-family symmetry groups which are compatible with explicit CP- 
violation. However we found another, quite remarkable feature in 3HDM: the presence of a 
Z 4 symmetry of Higgs-family transformations guarantees that the potential is explicitly CP- 
conserving. This is, of course, sufficient but not necessary condition for explicit CP-violation. 
Put in other words, explicit CP-violation is incompatible with the Higgs-family symmetry 
group Z 4 . 
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8.3 Further directions of research 

Certainly, our results do not provide answers to all symmetry-related questions which can be 
posed in 3HDM. Our paper should rather be regarded as the first step towards systematic 
exploration of all the possibilities offered by three Higgs doublets. Here are some further 
questions which deserve a closer study: 

• Continuous symmetry groups should also be included in the list. There exist only few 
Lie groups inside PSU(3): 17(1), £7(1) x £7(1), SU(2), SU{2) x £7(1), SO(3). The non- 
trivial question is which of these groups can be merged with some of the finite groups 
and with anti-unitary transformations (the case of abelian groups was analyzed in [2T]). 

• It is well-known that the vacuum state does not have to respect all the symmetries 
of the Lagrangian, so the finite symmetry groups described here can be broken upon 
electroweak symmetry breaking. What are the symmetry breaking patterns for each of 
these groups? Clearly, if the symmetry group is very small, then the vacuum state can 
either conserve it or break it, either completely or partially. But when the finite group 
becomes sufficiently large, there are two important changes: (1) the symmetry group 
cannot break down completely, as it would create too many degenerate vacua, which 
is not possible from the algebraic-geometric point of view; (2) some of the groups can 
never be conserved upon EWSB; the origin of this feature and some 3HDM examples 
were discussed in [30] . 

• What are possible symmetries of the potential beyond the unitary and antiunitary trans- 
formations? For example, the full reparametrization group of the 2HDM potential is 
GL(2,C) xi "LI rather than SU(2) x Z|, [20]. It means that a potential can be left in- 
variant by transformations which are neither unitary nor anti-unitary. Although these 
transformations played important role in the geometric constructions in the 2HDM orbit 
space, they did not produce new symmetry groups beyond what was already found from 
the unitary transformations. It would be interesting to check the situation in 3HDM. 
Unfortunately, the geometric method which worked well for 2HDM becomes much more 
intricate with more than two doublets, [TBI [11] . 

• It would also be interesting to see if the potential can have symmetries beyond repara- 
metrization transformations. In the case of 2HDM, this problem was analyzed in [3Tj. 
Although these additional symmetries cannot be extended to kinetic term, they could 
still provide useful information on the structure of the Higgs potential and properties of 
the physical Higgs bosons. 

In summary, we found all finite groups which can be realized as symmetry groups of Higgs- 
family or generalized-CP transformations in the three-Higgs-doublet model. Our list fl82|) is 
complete: trying to impose any other discrete symmetry group on the 3HDM Higgs potential 
will make it symmetric under a continuous group. 

This work was supported by the Belgian Fund F.R.S.-FNRS, and in part by grants RFBR 
11-02-00242-a, RFBR 10-01-0391, RF President grant for scientific schools NSc-3802.2012.2, 
and the Program of Department of Physics SC RAS and SB RAS "Studies of Higgs boson 
and exotic particles at LHC." 
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A Examples of the potentials with non-abelian Higgs- 
family symmetry group 

Here we collect the examples of potentials with a given symmetry group. These examples 
concern only non-abelian groups from the list (1521) because the abelian ones were already 
discussed in detail in [21]. We stress that these are just particular examples, and they of 
course do not exhaust all potentials symmetric under a given group. In fact, in the main text 
we described wider classes of such potentials than listed in this Appendix. 
Group Dq. Let us first denote the common part as 



Vx 



+ A13 



(4(h) + (0 2 2 )] - 774(0t 03 ) + An j (0 t 0i) 2 + (0 t 02) 2j + A33 ( t 03) 

(0100 + (4fo)} + A 12 (0 t i0 1 )(40 2 ) + Ai 3 \\4h\ 2 + 1403 



+ a; 



J3) 

|2 



where all coefficients are real and generic. Then, an example of the D 6 -symmetric potential is 



Vi 



D 6 



(0^0(0^0 - (0102)^02)1 + |A 3 |e^(0 t 1 3 )(0^03) + h.c. (84) 



where Ai is real and sin-0 3 ^ 0. This potential is symmetric under the group D G generated by 

(85) 














1 
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J 1 





I. -1 


1 












1 J 




V 








There are no other Higgs-family or generalized-CP transformations which leave this potential 
invariant. 

Group Dq x Z^. If in the previous example ( 1541) we set sin -03 = 0, then the potential 
becomes symmetric under D 6 x Zj. The extra generator, in addition to a and b, is just the 
complex conjugation. 

Group Dg x Zg. An example of the D§ x Zg-symmetric potential is 



D$xl 



V 1 + A 1 (0^0 1 )(0^0 2 ) + A 2 (0l0 2 ) 2 + h.c. 



(86) 



with the same V\ as before. Here both Ai and A2 are real (allowing them to be complex does 
not change the group, see the main text). The symmetry group is generated by 



If) 








°\ 






1 






—i 







1 





• 


V 





1 J 




^ I 








and the complex conjugation. 

Group A 4 x Zj. An example of the potential symmetric under A 4 x Z2 is 



V A 



4*U 



—m 
+A' 
+A 



>2<P2. 



+ 



+ A 



^l) + (4fc) + 



1 2 



(0100(0202) + (0 2 2 )(0303) + (0303)(0 t 10l)] + A" (\4<t>2? + |0 2 3 | 2 + \4<P 



(0l0 2 ) 2 + (0^0 3 ) 2 + (0^00 



+ h.c. 



(89) 
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with complex A. The symmetry group is generated by independent sign flips of the individual 
doublets, cyclic permutations of 0i, 2 , 03 and the exchange of any pair of doublet together 
with the simultaneous complex conjugation. An alternative form of this potential is 

n 2 



V, 



-m 



^0 



+ A 



400 + {4<h) + 



(90) 



+A' 



(0*00 (4^2) + (4h)(4h) + (0303) (0l0l] 



+A 



Rc 



(RC0I02) 2 + (Re0 2 3 ) 2 + (Re0 3 0i) 



+ A 



Im 



(WI02) 2 + (Im0 2 3 ) 2 + (Im0j0!) 



+A 



Rclm 



Re0{0 2 Im0j0 2 + Re0203 Im020 3 + Re030i 1UL030! 



Group S4 X Zg. An example of the potential symmetric under S4 x Zg is given by (1881) 
but with real A, or alternatively is given by (f90|) with AR e im = 0. The symmetry group is now 
generated by sign flips, all permutation of the three doublets, and the complex conjugation. 
Group (Z3 X Z3) x Z2 ~ A(54)/Z3. An example of potential with this symmetry is 

2 



V, 



A(54)/Z 3 



—m 

+4 

+A 2 
+A 3 



400 



^0 3 ) 



+ A c 



>l0i) + (0^02) 



403) 



(0I0!) 2 + (0 2 2 ) 2 + (0 3 3 ) 2 - (0100(0^02) " (0^0 2 )(0^03) ~ (4,<Pz)(4lK 



402 1 



+ 



4031 



+ 



40i I 



(0 t 1 2 )(0 t 1 3 ) + (0^0 3 )(0^0l) + (030l)(0302 



+ h.c. 



(91) 






with real m 2 , Ao, Ai, A2 and complex A3. All values here are generic. The symmetry group of 
this potential is (Z 3 x Z 3 ) xi Z 2 = A(54)/Z 3 , where A(54) is generated by 



(92) 



and Z 3 is the center of SU(3). 

Group (Z3 X Z3) xi (Z2 X Z^). The potential ( l9Ti) becomes symmetric under an antiunitary 
symmetry if 6 arg A 3 = 0. For example, if A 3 is real, the extra symmetry is complex conjugation. 

Group S(36) x Z^. The same potential (f9~Tj) becomes symmetric under the group £(36) x 
Z2 if the conditions (165]) are satisfied. An example of such potential is: 

2 



£(36) mZ^ 



m 2 (010! + 0202 + 0303 ) + A (4l4>l + 0202 + 0303 



+A 



402 



iir + |030i 



402 1 



(93) 



It is remarkable that this potential has only one "structural" free parameter A and that the 
term containing it reduces the full SU(3) symmetry group to a finite subgroup. 
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